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Abstract. A Goddard Earth Observing System Chemistry- 
Climate Model (GEOSCCM) simulation with strong tropi- 
cal non-orographic gravity wave drag (GWD) is compared 
to an otherwise identical simulation with near-zero tropical 
non-orographic GWD. The GEOSCCM generates a quasi- 
biennial oscillation (QBO) zonal wind signal in response to 
a tropical peak in GWD that resembles the zonal and clima- 
tological mean precipitation field. The modelled QBO has 
a frequency and amplitude that closely resembles observa- 
tions. As expected, the modelled QBO improves the simula- 
tion of tropical zonal winds and enhances tropical and sub- 
tropical stratospheric variability. Also, inclusion of the QBO 
slows the meridional overturning circulation, resulting in a 
generally older stratospheric mean age of air. Slowing of the 
overturning circulation, changes in stratospheric temperature 
and enhanced subtropical mixing all affect the annual mean 
distributions of ozone, methane and nitrous oxide. Further- 
more, the modelled QBO enhances polar stratospheric vari- 
ability in winter. Because tropical zonal winds are easterly in 
the simulation without a QBO, there is a relative increase in 
tropical zonal winds in the simulation with a QBO. Extrat- 
ropical differences between the simulations with and without 
a QBO thus reflect the westerly shift in tropical zonal winds: 
a relative strengthening of the polar stratospheric jet, polar 
stratospheric cooling and a weak reduction in Arctic lower 
stratospheric ozone. 


1 Introduction 

The quasi-biennial oscillation (QBO) is the leading mode 
of variability in the tropical lower and middle stratosphere 
(Baldwin et al., 2001). The QBO is characterized by a 
downward-propagating pattern of alternating easterly and 
westerly zonal winds in the equatorial region, with a pe- 
riod of approximately 28 months, and is driven by both grav- 
ity and planetary- scale waves (Lindzen and Holton, 1968; 
Holton and Lindzen, 1972). The zonal wind QBO induces 
changes in the tropical stratospheric circulation, affecting the 
concentrations of ozone and other trace constituents (Gray et 
al., 1989; Butchart et al., 2003; Tian et al., 2006). 

The phase of the QBO contributes to interannual vari- 
ability in the polar stratosphere. Holton and Tan (1980) and 
Lu et al. (2008) showed that the phase of the QBO modu- 
lates the strength of the Arctic vortex in mid-winter; the vor- 
tex is weakest during the easterly phase of the QBO. QBO 
phase-related differences in the strength of the polar vortices 
modulate polar ozone loss (Lait, 1989; Randel and Cobb, 
1994). Analogously to the Arctic response, the Antarctic vor- 
tex is generally weaker during the easterly phase of the QBO 
(Baldwin and Dunkerton, 1998). Also, the phase of the QBO 
modulates the polar stratospheric response to El Nino and La 
Nina events (Garfinkel and Hartmann, 2007; Hurwitz et al., 
2011a). However, as the QBO signal is intrinsic to the ob- 
servational record, the time-averaged impact of the easterly 
and westerly phases of the QBO on mean zonal wind, tem- 
perature and trace constituents is difficult to evaluate using 
atmospheric data. 
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A chemistry-climate model (CCM) is an ideal tool for un- 
derstanding the net impact of the QBO. Punge and Gior- 
getta (2008) quantified the net effect of the QBO on the late 
20th century stratospheric climate by comparing two CCM 
simulations: one without a QBO signal, and the other with 
zonal winds between 90hPa and lOhPa nudged to profiles 
taken at Singapore. The authors found that inclusion of this 
QBO signal impacted stratospheric zonal winds, temperature 
and ozone, mainly in the deep tropics. The vertical pattern of 
changes in tropical upwelling was consistent with changes 
in the distribution of trace species. However, the value of 
the conclusions reached by Punge and Giorgetta (2008) may 
be limited because the authors tested a “nudged” (i.e. pre- 
scribed) QBO. Nudging to observed zonal winds does not 
allow the tropical and mid-latitude stratosphere to adjust in- 
ternally to changes in, for example, the QBO phase, nor does 
it allow for full interaction between stratospheric ozone and 
climate. 

The Goddard Earth Observing System Chemistry-Climate 
Model (GEOSCCM) version 2 can be run with or with- 
out an internally generated QBO. In the formulation of 
the GEOSCCM evaluated by SPARC CCMVal (2010), both 
tropical lower stratospheric variability and the QBO ampli- 
tude were negligible. In contrast, a more recent model for- 
mulation (introduced by Hurwitz et al., 2011b) can internally 
generate a QBO with a realistic periodicity and amplitude, 
depending on the latitudinal structure of the non-orographic 
gravity wave drag (GWD). Comparing two simulations us- 
ing this formulation of the GEOSCCM, one with a QBO and 
another without, this paper quantifies the net effect of the 
modelled QBO on stratospheric climate and variability. Sec- 
tion 2 provides a brief description of the GEOSCCM, as well 
as the two above-mentioned simulations. The net effects of 
the modelled QBO on the mean and variance of stratospheric 
zonal winds, temperature, mean age, stratospheric circula- 
tion, ozone, nitrous oxide and methane are shown in Sect. 3. 
Section 4 provides a summary and brief discussion. 


2 Model and simulations 

This paper considers the net impact of the QBO in version 
2 of the GEOSCCM. The GEOSCCM couples the GEOS-5 
general circulation model (Molod et al., 2012) with a com- 
prehensive stratospheric chemistry module (Pawson et al., 
2008). The GEOSCCM performed well in the SPARC CCM- 
Val (2010) detailed evaluation of stratospheric processes. 
The present formulation of the GEOSCCM is the same as 
in Hurwitz et al. (2011b). In this formulation, the model has 
2° latitude x 2.5° longitude horizontal resolution and 72 ver- 
tical layers, with a model top at 0.01 hPa. Predicted distribu- 
tions of water vapour, ozone, methane, nitrous oxide, CFC- 
1 1, CFC-12 and HCFC-22 feed back to the radiative calcula- 
tions. 



Fig. 1. Background non-orographic gravity wave drag in the Q sim- 
ulation (shown as the solid black line) and N simulation (shown as 
the dashed black line). Reproduced with permission from Molod et 
al. (2012). 


In the GEOSCCM, tropical stratospheric zonal wind vari- 
ability depends on the details of the non-orographic GWD 
scheme. As the generation of non-orographic gravity waves 
often accompanies precipitation (e.g. convective and frontal 
systems; see Richter et al. (2010)), the latitudinal structure of 
the gravity wave spectrum is designed to mimic the structure 
of the climatological mean precipitation field (solid black 
line in Fig. 1; reproduced with permission from Molod et 
al., 2012). A 700 km wavelength is used for the tropical non- 
orographic waves to prevent an excessive downward propa- 
gation of the semi-annual oscillation into the lower strato- 
sphere, and thus contamination of the QBO signal. With 
a tropical peak in non-orographic gravity wave stress, the 
present model formulation generates a QBO signal in equa- 
torial zonal winds (see Sect. 3.1). Prior to this model formu- 
lation, the GEOSCCM did not generate a QBO: zonal winds 
in the equatorial lower stratosphere were generally easterly 
(see Fig. 2e and SPARC CCMVal, 2010). 

Two GEOSCCM simulations are used to assess the net 
stratospheric impacts of the QBO. The two simulations are 
identical, except for the magnitude of the non-orographic 
GWD stress applied in the deep tropics. The first simulation 
(hereafter “2”) has a strong peak in tropical GWD, which as 
described above, generates a QBO signal. The second simu- 
lation (hereafter “A”) has weak tropical GWD stress (dashed 
black line in Fig. 1), and thus does not have a QBO. 

Both simulations are 50 yr “time-slice” simulations with 
fixed climate forcings and annually repeating sea surface 
temperature (SST) and sea ice climatologies. Surface mix- 
ing ratios of the primary greenhouse gases and ozone- 
depleting substances are specified from 2005 concentrations. 
The SST and sea ice climatologies are composites of 10 El 
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Fig. 2. Equatorial (4° S-4° N) zonal wind (ms -1 ) in (a-b) MERRA reanalysis, 1979-2012, (c-d) the Q simulation, and (e-f) the N simu- 
lation. (a, c, e) lOyr time series, (b, d, f) frequency spectra for the entire span of the MERRA reanalysis and GEOSCCM simulations. Black 
dotted lines indicate the 30hPa level. 


Nino-Southern Oscillation (ENSO) -neutral years that span 
the satellite era (as in Hurwitz et al., 2011b). HadlSSTl SSTs 
and sea ice concentrations at 1° x 1° resolution (Rayner et 
al., 2003) are used to prepare the composites. Variability re- 
lated to the solar cycle and volcanic eruptions is not consid- 
ered. 

Modelled temperature and zonal wind fields are compared 
with the Modem Era Retrospective-Analysis for Research 
and Applications (MERRA). MERRA is a meteorological re- 
analysis, based on an extensive set of satellite observations 
and on the Goddard Earth Observing System Data Analysis 
System, version 5 (GEOS-5), from 1979 through the present 
(Rienecker et al., 201 1). The MERRA reanalysis has vertical 
coverage up to 0.1 hPa and, for this study, is interpolated to 
1.25° x 1.25° horizontal resolution. 

Additionally, simulated stratospheric mean age of air, 
ozone, nitrous oxide and methane are compared with obser- 
vational data sets. Simulated mean age of air is compared 
with profiles derived from CO2 and SF6 observations (An- 
drews et al., 2001; Engel et al., 2009). These observations 
were taken between 1986 and 2005. Simulated ozone (O3) 
and nitrous oxide (N2O) are compared with the 2004-2012 
climatology of the Aura Microwave Limb Sounder (MLS) 
version 3.3 O3 and N2O products (Froidevaux et al., 2008; 
Jiang et al., 2007; Livesey et al., 2008, 2011). Simulated 
methane (CH4) is compared with the UARS HALOE clima- 
tology for 1991-2002 (Russell et al., 1993; Grooss et al., 
2005). The UARS HALOE climatology is scaled by 1.02 
to reflect the increase in methane surface mixing ratio be- 
tween the mid-1990s and 2005 (see http://www.esrl.noaa. 
gov). More details about each data set are provided in Sect. 1 
of SPARC CCMVal (2010). 


3 Results 

3.1 Equatorial zonal winds 

The QBO in zonal winds is well simulated in Q. The left- 
hand panels of Fig. 2 show 10 yr time series of equatorial 
winds in the MERRA reanalysis (1990-1999; Fig. 2a), Q 
(Fig. 2c) and N (Fig. 2e). The right-hand panels of Fig. 2 
show the frequency spectra for equatorial winds for 1979- 
2012 in MERRA (Fig. 2b) and for the entire Q and N sim- 
ulations (Fig. 2d and f). The simulated QBO signal has re- 
alistic amplitude and periodicity (compare Fig. 2a and c). 
In Q , the simulated peak frequency is 27 months at 30hPa 
(Fig. 2d), with a secondary peak at 25 months, compared with 
28 months in the MERRA reanalysis (Fig. 2b). In the upper 
and middle stratosphere, the annual (12-month) and semi- 
annual (6-month) frequencies are well simulated. The annual 
frequency is weaker than observed in the lower stratosphere. 
Note that the simulated QBO signal is weaker than observed 
below 50 hPa (Hurwitz et al., 201 lb). 

Equatorial zonal winds are easterly throughout the N sim- 
ulation, without a QBO signal (Fig. 2e). That is, lower strato- 
spheric zonal wind variability is negligible. However, the an- 
nual and semi-annual frequencies are simulated in the middle 
and upper stratosphere (Fig. 2f). 

3.2 Annual mean impact of the QBO 

Inclusion of the QBO impacts stratospheric mean climate 
and variability. Figures 3, 5 and 6 show the 50 yr annual 
mean zonal wind, temperature, age of air, residual vertical 
and meridional velocities, methane, nitrous oxide and ozone 
fields in the N simulation (left-hand panels), Q simulation 
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(middle panels) and differences between Q and N (i.e. the 
net mean impact of the QBO) (right-hand panels). Fig- 
ure 4 shows QBO-related changes in stratospheric variabil- 
ity for each field. Annual and zonal mean zonal winds in 
the Q simulation are easterly in the tropical stratosphere 
and westerly in the extratropics (Fig. 3b), in good agree- 
ment with the MERRA reanalysis (black contours). In N , 
tropical stratospheric easterlies are on average stronger than 
in the MERRA reanalysis (Fig. 3a); thus, inclusion of the 
QBO represents both a relative zonal wind increase in the 
tropics (10-20 ms -1 ) and an improvement in the simulated 
mean comparison with MERRA (compare Fig. 3a and b). As 
expected, zonal wind variability increases by 3-5 times in 
2, as compared with N , throughout the tropical stratosphere 
(Fig. 4a). 

The modelled stratospheric temperatures generally agree 
well with MERRA (Fig. 3d and e). A modest cold bias 
around 60° S may reflect the overly strong polar jet and de- 
layed breakup of the Antarctic vortex in austral spring (Hur- 
witz et al., 2010). Inclusion of the QBO warms the tropi- 
cal stratosphere and Arctic lower stratosphere, but cools the 
high-latitude upper stratosphere by ~ 1 K (Fig. 3f). These 
temperature changes reflect a slowing of the stratospheric 
overturning circulation, as discussed below. QBO-related in- 
creases in temperature variability maximize in the middle 
stratosphere, in the deep tropics, with an additional lobe of 
increased variability around 30° S (Fig. 4b). 

Inclusion of the QBO acts to slow the stratospheric circu- 
lation. Age of air is an indicator of the strength and struc- 
ture of the meridional overturning (i.e. Brewer-Dobson) cir- 
culation (SPARC CCMVal, 2010). In the GEOSCCM, mean 
age is an inert tracer that measures the time since a parcel 
of air has left the troposphere. Because most air enters the 
stratosphere in the tropics, mean age values are lowest in the 
tropical lower stratosphere and highest in the high-latitude 
upper and middle stratosphere (Fig. 5a and b). Q — N differ- 
ences in annual mean age are generally positive throughout 
the stratosphere i.e., the QBO slows the overturning circula- 
tion (Fig. 5c). The largest increases in mean age (~ 0.3 yr) 
are seen in the tropical stratosphere around lOhPa. Figure 7 
shows observed and simulated mean age-of-air profiles in the 
deep tropics. While the mean age in N falls within the ob- 
servational error, the older tropical mean age in the Q sim- 
ulation is in better agreement with CO 2 and SF6 observa- 
tions. QBO-related changes in mean age variability, similarly 
to temperature variability, peak in the Southern Hemisphere 
tropical middle stratosphere (Fig. 4c). Decreased mean age 
in the lower stratosphere is consistent with slowing of the 
overturning circulation: less older air is advected downward 
into the mid- latitudes. 

QBO-related changes in the residual circulation are con- 
sistent with the changes in mean age of air. In the upper 
stratosphere and polar regions, the net impact of the QBO 
is dominated by the slowing of the meridional circulation. At 
and above lOhPa, Q — N differences in (residual vertical 


velocity) oppose the climatology: negative differences in 
the tropics (weakened upwelling), and positive differences 
in the extratropics (weakened downwelling; compare Fig. 5e 
and f). In Fig. 5g-i, positive 1 ;* (residual meridional veloc- 
ity) values indicate poleward motion in both hemispheres. 
Negative differences in 1 ;* (e.g. around 3hPa) in the trop- 
ics and in the uppermost polar stratosphere reflect weakened 
poleward transport (Fig. 5i). In the lower and middle strato- 
sphere, the net impact of the QBO is to enhance poleward 
transport. Positive between 10 and 70hPa, in the tropics 
and mid-latitudes, is consistent with enhanced mixing in this 
region. Note that the 1 ;* differences are largest in the South- 
ern Hemisphere and coincident with the peak change in 1 ;* 
variability (Fig. 4e). Below 70hPa, tropical 1 ;* changes are 
positive. This set of * changes is consistent with the nar- 
rowing but strengthening of upward transport in the tropical 
lower stratosphere (Fig. 5f). 

Both slowing of the stratospheric overturning circulation 
and enhanced lower stratospheric mixing affect the distribu- 
tion of methane. Since methane is mainly transported into 
the stratosphere in the deep tropics, the distribution of strato- 
spheric methane is inverted as compared with mean age: 
the highest concentrations are found in the tropical lower 
stratosphere, while the lowest values are found in the po- 
lar upper stratosphere. The structure of the simulated trop- 
ical and subtropical methane distribution matches the scaled 
UARS HALOE climatology, but with a high bias (Fig. 6a 
and b). This bias results from the stronger than observed 
transport in the GEOSCCM (SPARC CCMVal, 2010). Inclu- 
sion of the QBO decreases tropical and upper stratospheric 
methane mixing ratios, reflecting slowing of the stratospheric 
circulation (Fig. 6c). Given the negative latitudinal gradient 
in methane concentrations (Fig. 6b), increased lower strato- 
spheric methane suggests enhanced meridional mixing (i.e. 
weakening of the subtropical transport barrier); this effect is 
further discussed below. Inclusion of the QBO generally in- 
creases methane variability, with the largest increases in the 
tropics and subtropical Southern Hemisphere (Fig. 4f). 

QBO-related changes in N 2 O provide further evidence 
for enhanced subtropical mixing. Similarly to methane, the 
highest N 2 O concentrations are found in the tropical lower 
stratosphere, while the lowest values are found in the upper 
stratosphere (Fig. 6d and e). The simulated N 2 O distribution 
matches the MLS climatology, but with a high bias in the 
tropical lower stratosphere, similar to the methane bias. In- 
clusion of the QBO weakens the meridional gradient of N 2 O: 
decreasing equatorial N 2 O (exceeding 30ppbv) and increas- 
ing N 2 O at subtropical latitudes (up to lOppbv) (Fig. 6f). 
Figure 8 shows the annual cycle of Q — N differences in the 
subtropical N 2 O gradient. Mean N 2 O differences between 
10 and 40° latitude serve as a proxy for the strength of the 
subtropical gradient (i.e. the strength of the subtropical mix- 
ing barrier). Negative differences indicate a weakening of 
the subtropical N 2 O gradient, due to enhanced mixing (Dou- 
glass et al., 1999). The largest changes in subtropical N 2 O 
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Fig. 3. (a-c) Zonal wind (ms -1 ) and (d-f) temperature (K), for (a, d) N simulation zonal and annual mean and (b, e) Q simulation zonal 
and annual mean; black contours indicate MERRA reanalysis mean for 2000-2010. (c, f) Q — N mean differences; white contours indicate 
zero difference; black Xs indicate changes significant at the 95 % level in a two-trailed t test. 
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Fig. 4. Q/N changes in variance. White contours indicate where the variance ratio is equal to one. Black Xs indicate changes significant at 
the 95 % level in an / test. 


are centred around lOhPa (consistent with the increases in 
subtropical u* and methane), in winter months. The QBO 
has a larger impact on subtropical N 2 O in the SH (Fig. 8a) 
than in the NH (Fig. 8b), consistent with the relatively larger 
enhancements in variability in the SH (Fig. 4). 


The net impact of the QBO on ozone reflects warm- 
ing of the tropical stratosphere, enhanced subtropical mix- 
ing and slowing of the stratospheric overturning circula- 
tion. The annual mean ozone mixing ratio maximizes at 
approximately lOppmv in the deep tropics, at lOhPa, and 
decreases with latitude (Fig. 6g and h). Simulated ozone 
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Fig. 5. As for Fig. 3, but for (a-c) mean age (years), (d-f) w*(ms 1 ) and (g-i) w(ms 1 ). Note that values of in the Southern Hemisphere 

are multiplied by —1, such that positive values indicate poleward motion in both hemispheres. 
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Fig. 6. As for Fig. 3, but for (a-c) methane (ppmv), (d-f) nitrous oxide (ppbv) and (g-i) ozone (ppmv). Black contours in (a) and (b) indicate 
the scaled HALOE climatology; black contours in (d-e) and (g-h) indicate MLS climatologies of nitrous oxide and ozone. 


in the Q simulation is generally in agreement with the 
MLS ozone climatology (Lig. 6h). Thus, inclusion of the 
QBO improves model performance; the peak in ozone mix- 
ing ratio around lOhPa is ~ 0.8 ppmv too high in the 


N simulation. Relative warming of the tropical middle 
stratosphere in Q contributes to the negative ozone differ- 
ences in the deep tropics (Ligs. 3f and 6i). However, at 
the Equator at lOhPa, the ratio between temperature and 
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Fig. 7. Equatorial (10° S-10° N) mean age-of-air (years) profiles in 
the N simulation (red), Q simulation (blue), and from CO 2 (black 
diamonds) and SF^ (black squares) observations. 
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Fig. 8. Q — N differences in the subtropical N 2 O gradient (ppbv), as 
a function of month and altitude in (a) the Southern Hemisphere and 
(b) the Northern Hemisphere. Negative differences indicate weak- 
ening of the subtropical mixing barrier (i.e. enhanced subtropical 
mixing). White contours indicate zero differences. 


ozone differences (i.e. AT/ AO 3 = 0.72K/0.8ppmv = 0.9) 
does not agree with the ratio of ~ 6.7 as determined by Oman 
et al. (2010). This result suggests that other mechanisms con- 
tribute to the change in tropical ozone. Decreased tropical 
ozone at lOhPa is consistent with both (1) strengthening of 
vertical upwelling below lOhPa (Fig. 5f) and (2) enhanced 
mixing with mid-latitude air with relatively lower ozone con- 
centrations (Fig. 6h and i). 

3.3 High-latitude impact of the QBO in winter 

In certain months, in the Q simulation, equatorial zonal 
winds are positively correlated with zonal winds in the po- 
lar stratosphere. That is, the GEOSCCM reproduces the 
Holton- Tan (1980) relation between the phase of the QBO 
and polar vortex strength: in winter, the polar vortices are 
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Fig. 9. 50 hPa zonal wind correlation with equatorial zonal winds, 
as a function of month and latitude, in (a) MERRA and (b) the Q 
simulation. White contours indicate zero correlation. 


relatively stronger during the westerly phase of the QBO 
in the MERRA reanalysis (Fig. 9a). Figure 9b shows that, 
in the GEOSCCM, equatorial zonal wind correlations with 
zonal winds at 60° S are strong in, for example, September 
(Fig. 10), and at 60° N in January (Fig. 12). 

The net impact of the QBO on zonal winds and tempera- 
ture (i.e. Q — N ) mimics the differences between the westerly 
and easterly phases of the QBO (i.e. in the Q simulation), 
at Southern Hemisphere high latitudes in September. Inclu- 
sion of the QBO strengthens the Antarctic jet (around 60° S) 
by 5-10 ms -1 (Fig. 10a). An equal strengthening is seen in 
Fig. 10b, which shows the zonal wind differences between 
the five Septembers in Q with the highest (i.e. QBO-W) 
equatorial zonal wind values at 30hPa and the five Septem- 
bers in Q with the lowest (i.e. QBO-E) equatorial zonal wind 
values at 30 hPa. Note that the QBO-W — QBO-E differences 
are generally not statistically significant in the extratropics, 
due to the small size of the QBO-W and QBO-E composites. 

Extratropical stratospheric temperatures decrease in 
September, consistent with the stronger zonal winds (i.e. 
by the thermal wind relation). Temperature differences are 
negative poleward of 45° S, both for Q — N and QBO-W 
— QBO-E, with peak differences of 5 K in the Antarctic 
(Fig. 10c and d). Tropical and mid-latitude Q — N differences 
reflect changes in stratospheric circulation and mixing (see 
Sect. 3.2). 

Strengthening and cooling of the Antarctic polar vortex in 
September reflects the relative shift toward westerly tropical 
zonal winds, shifting the zero-wind line equatorward and de- 
creasing planetary wave breaking at high latitudes. The mag- 
nitude of eddy heat flux at 40-80° S and lOOhPa diagnoses 
tropospheric planetary wave driving (Newman et al., 2001). 
The red line in Fig. 11 shows that eddy heat flux is reduced 
in the Q simulation relative to the N simulation, in August 
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Fig. 10. September (a-b) zonal wind (ms -1 ), (c-d) temperature (K) and (e-f) ozone (ppmv) differences in the Southern Hemisphere, (a, c, 
e) Q — N mean differences; (b, d, f) QBO-westerly — QBO-easterly differences in the Q simulation. White contours indicate zero difference. 
Black Xs indicate differences significant at the 95 % level in a two-trailed t test. 



Fig. 11. Q — N differences in eddy heat flux magnitude at 40°-80° 
latitude, 100 hPa in the Northern Hemisphere (blue line) and South- 
ern Hemisphere (red line). Negative differences indicate a relative 
reduction in tropospheric planetary wave driving. The case shown 
in Fig. 10 is highlighted by the red circle. The case shown in Fig. 12 
is highlighted by the blue circle. 

through October. Note that the annual mean eddy heat flux 
difference is negligible. 

Equivalent zonal wind and temperature differences are 
seen in the Northern Hemisphere extratropics in January. 
Q — N differences reflect a relative shift toward the west- 
erly phase of the QBO: a relative strengthening of the Arctic 
stratospheric jet (Fig. 12a and b) and cooling throughout the 


extratropical stratosphere (Fig. 12c and d). In this case, zonal 
wind and temperature differences are larger in response to 
QBO phase as compared with the net impact of the QBO. 
Strengthening and cooling of the Arctic vortex is caused by 
a dramatic reduction in planetary wave driving in January 
(Fig. 11, blue dot). Negative ozone differences in the Arc- 
tic lower stratosphere (Fig. 12e and f), though weak, hint at 
enhanced chemical ozone depletion, due to the cooler and 
more isolated polar air mass in January. The small size of the 
QBO-W and QBO-E composites likely weakens the statisti- 
cal significance of the ozone differences. 

4 Discussion 

Inclusion of the QBO in a chemistry-climate model makes a 
significant difference to the mean stratospheric climate and 
variability. As expected, the addition of a QBO significantly 
enhances tropical variability. Extratropical zonal wind and 
temperature variability in winter is also enhanced. While the 
QBO is by definition an oscillating phenomenon, the multi- 
decadal mean of the modelled QBO modifies the mean cli- 
mate. In particular, the following can be concluded: 
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Fig. 12. January (a-b) zonal wind (ms -1 ), (c-d) temperature (K) and (e-f) ozone (ppmv) differences in the Northern Hemisphere, (a, c, e) 
Q — N mean differences; (b, d, f) QBO-westerly — QBO-easterly differences in the Q simulation. White contours indicate zero difference. 
Black Xs indicate differences significant at the 95 % level in a two-trailed t test. 


1. Addition of a QBO signal affects the mean strato- 
spheric climate. In the GEOSCCM, inclusion of a 
QBO signal slows the meridional overturning circu- 
lation, leading to older mean age of air throughout 
the stratosphere and affecting the distribution of trace 
species such as ozone, methane and nitrous oxide. 
The subtropical mixing barrier is weakened, enhanc- 
ing methane and nitrous oxide in the extratropical 
lower stratosphere, and contributing to the reduced 
peak ozone mixing ratio at 1 0 hPa. The polar vortices 
are strengthened, particularly in winter. The net dy- 
namical and transport impacts of the QBO generally 
improve model performance. 

2. The net impact of the QBO depends on both the base- 
line zonal wind field and the relative change in tropical 
zonal winds. In the GEOSCCM, the baseline tropical 
zonal winds are easterly. In the simulation with a QBO 
signal, there is thus a relative increase in zonal winds 
in the tropical lower and middle stratosphere. Extrat- 
ropical differences between simulations with and with- 
out a QBO reflect this shift toward the westerly zonal 
winds in the tropical stratosphere: a relative cooling 
and strengthening of the polar vortices, and a weak re- 


duction in Arctic lower stratospheric ozone. The an- 
nual mean impact of the QBO on the polar strato- 
sphere is larger in the GEOSCCM (up to 12ms -1 , in 
both hemispheres) than in the MAECHAM4-CHEM 
CCM (no significant zonal wind differences) (Punge 
and Giorgetta, 2008), likely reflecting larger tropical 
zonal wind differences between the QBO and “no 
QBO” GEOSCCM simulations and/or increased sta- 
tistical robustness due to the relatively greater length 
of the GEOSCCM simulations. 


The QBO has a robust net impact on the mean strato- 
spheric climate and trace gas distributions. In the case of the 
GEOSCCM, inclusion of the QBO yields better agreement 
between the simulated fields and climatological averages de- 
rived from a meteorological reanalysis and satellite data sets 
of ozone, methane and nitrous oxide. While it is difficult to 
generate a QBO signal internally in a global climate model, a 
model with the ability to simulate a QBO presents significant 
advantages in predicting the future evolution of the strato- 
sphere: climate change is likely to modify the GWD, which 
will remotely modify stratospheric climate and variability. 
Capturing this feedback will require a modelled QBO forced 
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by an interactive GWD source spectrum, driven by changes 
in, for example, convective precipitation. 
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